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ABSTRACT 

In this spectroscopic study of infant massive star clusters, we find that continuum emission from 
ionized gas rivals the stellar luminosity at optical wavelengths. In addition, we find that nebular 
line emission is significant in many commonly used broad-band Hubble Space Telescope (HST) filters 
including the F814W /-band, the F555W U-band and the F435W B-band. Two young massive 
clusters (YMCs) i n the nearby starburst galaxy NGC 4449 were targeted for follow-up spectroscopic 



observations after Reines et al. (2008a) discovered an F814W /-band excess in their photometric 
study of radio-detected clusters in the galaxy. The spectra were obtained with the Dual Imaging 
Spectrograph (DIS) on the 3.5 m Apache Point Observatory (APO) telescope]^] and have a spectral 
range of ~3800-9800 A. We supplement these data with HST and Sloan Digital Sky Survey (SDSS) 
photometry of the clusters. By comparing our data to the Starburst99 and GALEV evolutionary 
synthesis models, we find that nebular continuum emission competes with the stellar light in our 
observations and that the relative contribution from the nebular continuum is largest in the U- and 
/-bands, where the Balmer (3646 A) and Paschen jumps (8207 A) are located. The spectra also 
exhibit strong line emission including the [S III] AA9069,9532 lines in the HST F814W /-band. We 
find that the combination of nebular continuum and line emission can account for the F814W /-band 



excess previously found by Reines et al.| (|2008aJ) . In an effort to provide a benchmark for estimating 
the impact of ionized gas emission on photometric observations of young massive stellar populations, 
we compute the relative contributions of the stellar continuum, nebular continuum, and emission lines 
to the total observed flux of a 3 Myr-old cluster through various HST filter/instrument combinations, 
including filters in the Wide Field Camera 3 (WFC3). We urge caution when comparing observations 
of YMCs to evolutionary synthesis models since nebular continuum and line emission can have a 
large impact on magnitudes and colors of young (< 5 Myr) clusters, significantly affecting inferred 
properties such as ages, masses and extinctions. 

Subject headings: galaxies: individual (NGC 4449) — galaxies: starburst — galaxies: star clusters — 
H II regions — ISM: lines and bands 



1. INTRODUCTION 

Massive star clusters are an important mode of star 
formation, having an impact on a wide range of galaxy 
properties. However, the earliest stages of these clusters 
are notoriously challenging to study since the youngest 
clusters are still enshrouded in remnants of their gaseous 
and dusty birth cocoons. Properly accounting for the ef- 
fects of gas and dust on observations of young massive 
clusters (YMCs) is nontrivial, rendering it difficult to un- 
derstand their formation and earliest phases of evolution. 
Some of these difficulties, such as interstellar extinction, 
are commonly dealt with in studies of young star clus- 
ters. While it is often not possible to completely dis- 
entangle the effects of extinction, t he general effects are 



well-known and accounted for (e.g. Calzetti e t al. 1994 
Whitmore et al.||1999| |Johnson et al.||1999| |Reines et al 
2008b! - 



On the other hand, the treatment of ionized gas emis- 
sion can be quite complicated and the effects of nebu- 
lar emission on broad-band observations of YMCs are 
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not always obvious. Emission lines detected via spec- 
troscopy or narrow-band imaging, however, have proved 
useful for both identifying young (< 10 Myr) clusters 
(e.g. Ha) and serving a s valuable diagnost i cs on their 
physical conditions (e.g. |Eissner et al. || 1969 | |Castan eda 
et al.||1992| IKewley & Dopita||2002p. if these lines 



are 



strong enough, as is typically the case for massive star 
forming regions, the composite emission can even affect 
the integrated broad-band photometry (e. g. Johnson et 



aL] |1999l |Anders fc Fritze-v. Alvensleben|[2003p . How 
ever, nebular emission also takes on a more subtle form 
via continuum emission (e.g. free-free, free-bound). In 
fact, discontinuities or "jumps" in the nebular contin- 
uum (e.g. the Balmer and Paschen jumps) can be strong 
enough to serve a s diagnostics in their own right (see 



Guseva et al.|2006 and references therein) . This nebular 



continuum emission can also severely affect t he broad 
band photometry of very young clu sters (e.g. Leitherer 
|fc Heckman1p95| |Molla et al.|[2009 1. 

In this paper we investigate the impact of nebular con- 
tinuum and line emission on observations of very young 
massive star clusters. We carry out a detailed multi- 
wavelength study of two YMCs that are still partially 
embedded in their birth material. The clusters were 
targete d for spectroscopic follow-up after |Reines et al. 
(2008a) found an 7/STF814W /-band excess in their pho- 
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Fig. 1. — An SDSS 3-color image of NGC 4449 (RG B=uqi) using a logar ithmic scaling. The two YMCs for which we have spectra are 
labeled (Sources 26 and 30 according to the notation of|Reines et al.|2008a|. The field of view is approximately 3' X 3'. 



tometric study of radio-detected c lusters in the nearby 
dwarf starburst galaxy NGC 4449. |Reines et aE] ( |2008a I 
considered several possible explanations for this 1 -band 
excess, including red supergiants, emission lines, ther- 
mal emission from hot dust, continuum emission from a 
sub-population of dee ply embedded stars , and Extended 
Red Emission (ERE; [Witt fc Vijh|[2004| . With the ob- 
servations presented here, we have determined that the 
combination of nebular continuum and lin e emission can 
account for the /-band excess found by |Reines et al.| 
(2008a I. Consequently, the chief goal of this study is 



to gain a better empirical understanding of the effects of 
nebular emission (both continuum and line) associated 
with extremely young massive star clusters. 

This paper is organized as follows: The data are de- 
scribed in ^J2j In Sj3]we estimate the physical properties of 
the clusters from their Ha emission. We model the spec- 
tral energy distributions (SEDs) in Sj4] and investigate 
the impact of nebular emission (line and continuum) on 
broad-band photometry in ^j5] A summary of our main 



conclusions is given in Sj6] 

2. DATA 

A three-color SDSS image of NGC 4449 is shown in 
Figure [T] with our target Y MCs labeled. Sourc e 26 (ac- 
cording to the notation of |Reines et al. 2008a) was se- 
lected for its brightness (to get a high signal-to-noise 
spectrum) and Source 30 was targeted for its large I- 
band excess - the lar gest of all of the clusters studied in 
Reines et al. ( 2008a|. We adop t a dist ance of 3.9 Mpc 
to 1NGC 4449 as in|Reines et al.| (f2008a[ ) , consistent with 
the results of Annibali et al. 1 2008p . At this distance, 
1" - 19 pc. 

2.1. Spectroscopy 

Spectra of the two YMCs were obtained using the 3.5 m 
APO telescope. Source 26 (RA = 12 h 28 m 13.86 s , DEC = 
+44° 07' 10.4") was observed on the night of 2008 April 
13 under clear conditions and Source 30 (RA = 12 h 28 m 
16.02 s , DEC = +44° 06' 29.3") was observed on 2009 
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Fig. 2. — APO/DIS flux calibrated spectra of our target YMCs. HST and SDSS broad-band photometry are plotted as dots and the 
approximate widths of the filters are shown as gray bars. The spectra exhibit many strong emission lines as well as nebular continuum 
emission as is evident by the Paschen jump at 8207 A. The "W-R bump" at ~ 4650 A signifying the presence of Wolf-Rayet stars is also 
clearly visible in the spectrum of Source 26. 



May 24 through intermittent clouds. The observations 
were made using the red and blue channels of the Dual 
Imaging Spectrograph (DIS) in the low-resolution mode 
with the B400 and R300 gratings. The linear dispersions 
for the blue and red channels are 1.83 and 2.31 A pix -1 , 
respectively, and the central wavelengths are —4400 and 
-7500 A. We used the l'/5 x 360" slit oriented along the 
parallactic angle and centered on the YMCs. The result- 
ing spectra have a wavelength range of -3800-9800 A 
and a resolution of —7 A. The total exposure time on 
Source 26 was 30 minutes (2 x 15 minutes), resulting 
in signal-to-noise ratios of —60 and —45 in the blue and 
red continua, respectively. Source 30 was observed for 
a total of 2.5 hours (5 x 30 minutes) and the signal- 
to-noise ratio of both the blue and red continua is —5. 
Multiple exposures of each cluster were taken to allow 
for the removal of cosmic rays. The white d warf primary 



spect rophot omet ric standard star GD 153 ( |Bohlin et al 
1995 1 was observed for flux calibration and a He-JNe-Ar 



arc lamp spectrum was taken for wavelength calibration. 

The data were reduced using a combination of IRAF 
and custom IDL routines. The two-dimensional images 
were processed using standard procedur es in IRAF in - 
cluding bias subtraction and flat-fielding (Ma ssey |1997 |. 
IRAF was also used to extract the one-dimensional spec- 
tra, subtrac t the background, an d apply the wavelength 



calibration (Massey et al. 19921. The spectra were ex- 



tracted from an aperture window of 1 1 pixels (4. 4) in the 
spatial dimension. The remaining data reduction and 
analysis was performed in IDL. Flux calibration was ob- 



tained by comparing the observed spectrum of GD 153 to 
the model spectrum obtained from the Space Telescope 
Science Institute's Calibration Data Base System. 2 At- 
mospheric extinction was corrected for using the stan- 
dard extinction curve for APO. Finally, the blue and red 
spectral regions were combined and put on a common 
wavelength scale of 2 A pix -1 . The blue and red spec- 
tra were averaged in an overlap region approximately 50 
pixels wide and centered at A — 5370A. 

After a comparison with the photometry (see below), 
the flux calibrated spectra were multiplied by scale fac- 
tors of 1.94 (Source 26) and 1.84 (Source 30) to account 
for loss of light in the slit. We scaled the spectra to match 
the HST F550M (narrow U-band) data since this filter 
has a negligible contribution from emission lines. Fig- 
ure [2] shows the final reduced spectra of Sources 26 and 
30. The spectra exhibit many strong emission lines and 
nebular continuum emission as is evident by the Paschen 
jump at 8207 A. The Wolf-Rayet blended emissio n fea- 
ture at — 4650 A known as the "W-R bump" (e.g. |Conti| 
1991 Schaerer fc Vacca|1998 1 is also clearly visible in the 
spectrum of Source 26. 

2.2. Photometry 

HST observations of NGC 4449 taken with the Ad- 
vanced Camera for Surveys (ACS) and the Wide Field 
and Planetary Camera 2 (WFPC2) were obtained from 
the archive. These data are described in detail in IReinesl 

2 http:/ /www. stsci.edu/hst/observatory/cdbs/calspec. html 
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TABLE 1 

Archival HST Observations of NGC 4449 



Filter 


Instrument 


Description 


Exp. Time (s) 


Proposal ID 


PI 


F170W 


WFPC2 


UV 


400 X 2 


6716 


T. Stecher 


F336W 


WFFC2 


WKFC2 U 


2bU x 2 


bVlb 


1 . Stecner 


F43bW 


ACS/WKC 


jonnson b 


3bbU 


!Ub8b 


A. Aloisi 


FbbOM 


ACS/WKC 


IN arrow V 


120U 


!Ub22 


D. uaizetti 


Fbb8iN 


ACS/WKC 


Ho + fNlllAb'b84 


lb39 


10b8b 


A. Aloisi 


FbbON 


ACS/WKC 


[Nlf|Abb84 


18b0 


10b22 


u. uaizetti 


F814W 


ACS/WKC 


LSroad 1 


2060 


10b8b 


A. Aloisi 




TABLE 2 
Flux Densities of the YMCs 



Filter A (A) Source 26 Source 30 



HST 



F170W 


1821 


1.3(0.1) x 10- 


14 


2.2(0.7) x 1(T 


16 


F336W 


33b9 


4.1(0.4) x 10~ 


15 


1.4(0.1) x 10" 


16 


F43bW 


4317 


2.2(0.2) x 10" 


IB 


7.8(0.8) x 10- 


17 


FbbOM 


bb81 


9.2(0.9) x 10- 


16 


3.0(0.3) x 10- 


17 


F6b8N 


6b84 


1.6(0.2) x 10" 


14 


1.1(0.1) x 10- 


15 


F660N 


6b99 


3.4(0.3) x 10" 


15 


2.6(0.3) x 10" 


16 


F814W 


8061 


5.4(0.b) x 10- 


16 


3.0(0.3) x 10- 


17 



SDSS 



(( 


3b43 


4.3(0.4) X 10 


-15 


1.2(0.1) X 10- 16 


g 


4770 


3.0(0.3) x 10' 


-15 


1.6(0.2) x 10" 16 


r 


6231 


1.8(0.2) x 10' 


-15 


8.3(0.9) x 10~ 17 


i 


762b 


4.9(0.b) x 10" 


-16 


1.8(0.2) x 10" 17 


z 


9134 


4.0(0.4) x 10 


-16 


1.7(0.3) x 10- 17 


Note. - 


- Units of flux density are 


erg s" 


1 cm - 2 A -1 . Uncer- 


tainties are 


shown in 


parentheses. 








Fig. 3. — An HST 3-color image of NGC 4449 in the vicinity of 
Sources 26 and 30. Red, green and blue correspond to the F6b8N 
(Ha), F814W (broad /), and FbbOM (narrow V) filters. The circles 
indicate the phot ometric apertures used on the HST and SDSS 
imaging (see §2jh. 



et al. (l2008at) and summarized in Table [I] Here we 
use WFPC2 F170W and F336W images, as well as 
ACS Wide Field Channel (WFC) images through the 
F435W, F550M, F658N, F660N and F814W filters. The 
broad-band filters (F170W, F336W, F435W, F550M and 
F814W) span ultraviolet through near-IR wavelengths. 
The narrow-band filters, F658N and F660N, capture 
Ha+ [N II] A6584 and [N II] A6584 emission, respectively. 
We have also obtained ugri z imaging of NGC 4449 from 



the SDSS Data Release 6 (York et al. 2000 Adelman 



McCarthy et al.||2008[ ). A three-color (R<J,B=ugi) image 



of the galaxy is shown in Figure [T] with our target YMCs 
labeled. 

Photometry of the YMCs was performed using SUR- 



PHOT ( Reines et aL]|2008a| ), a custom IDL program al- 
lowing tor consistent apertures and background annuli 
across multiple wavebands. Circular apertures of radii 
3'.'3 and l'/7 were used for Sources 26 and 30, respectively. 
Background levels were determined in annuli with inner 
and outer radii equal to 1.75x and 2.5x the aperture 
radii. Aperture sizes were chosen to include essentially 
all of the light from the clusters in all bands, and were 
mostly influenced by the (sccing-limited) resolution of 
the ground-based SDSS images and the intrinsic cluster 
sizes. Figure [3] shows a three-color HST image of the 
clusters with the photometric apertures overlaid. Flux 
densities for the YMCs are given in Table [2] and plotted 
in Figure [2] along with the spectra. 

2.3. A Comparison of the Spectroscopic and 
Photometric Data 

2.3.1. Broad-band Flux Densities 

In order to compare our spectroscopic observations 
with the broad-band imaging data, we simulate "DIS 
photometry" by convolving the YMC spectra with the 
total system throughput curves for the g, F550M, r, i, 
and F814W filters. 3 Figure [5] shows the broad-band flux 
densities obtained from the imaging data (listed in Ta- 
ble [2]), as well as the flux densities obtained from the 

3 The observed spectra do not cover the entire wavelength ranges 
of the F43bW and z filters. 
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Fig. 4. — Broad-band flux densities obtained from the imaging 
data (diamonds) and from convolving the observed spectra with the 
total system throughput curves for the filters shown (plus signs). 
The sizes of the plot symbols reflect the errors. The spectroscopic 
data has been scaled to match the F550M photometry to account 
for loss of light in the slit. The discrepancy between the g flux 
densities for Source 30 are most likely due to an incorr ect absolute 
flux level in the blue part of the spectrum (see j ]2.3.f \ . 

spectroscopy. 

Figure [3] illustrates that the flux densities obtained 
from convolving the spectra with the filter throughput 
curves match the the photometry quite well, especially 
for Source 26 which has a high signal-to-noise spectrum. 
However, the SDSS g flux densities for Source 30 are 
clearly discrepant, with the spectroscopic value signifi- 
cantly higher than the photometric value. We believe 
this is due to an incorrect absolute flux level in the blue 
part of the spectrum for two reasons. First, we expect the 
F435W photometric data point to lie above the contin- 
uum in Figure [2] since this filter contains strong emission 
lines (e.g. HTand US). This effect can be seen for Source 
26 in Figure [2] (i.e. the F435W data point lies above the 
continuum in the spectrum), and the emission lines are 
even stronger relative to the continuum for Source 30. 
The second reason we think the blue end of the spec- 
trum of Source 30 has an incorrect flux level is that no 
single model SED can reproduce the shape of the contin- 
uum in the spectrum. There is no combination of age and 
extinction that can simultaneously produce such a steep 
blue spectrum and flattened red spectrum. Also, we show 
later in f|4]that the best-fitting SED to the broad-band 
photometric data has a much flatter slope than the blue 
part of the spectrum of Source 30. 

It is not absolutely clear why the spectrum of Source 
30 has an incorrect absolute flux level in the blue, al- 
though it is most likely due to the bad weather (clouds) 
during the observations. We note that besides the g- 
band, the flux densities derived from the spectroscopy 
are consistent with the HST and SDSS photometry, and 
that the error in the absolute flux level in the blue part 
of the spectrum of Source 30 does not affect our analysis 
or conclusions. 

2.3.2. Ha Emission 

In studies of extragalactic massive star forming regions, 
narrow-band imaging is commonly used to estimate the 
emission of Ha, which in turn is used to infer physi- 



cal properties such as ionizing luminosity, age and mass. 
Since we have spectroscopy and narrow-band imaging of 
our target clusters, we compare the flux and equivalent 
width of Ha from both data sets to see how well the 
measurements agree. 

Using the HST photometry, the Ha flux, Fn a , and 
equivalent width, Wn a , arc given by 



/F658N — /cont ) AAf65x\ 



^/F660N — /iontP) AAp660N 



(1) 



and 



F, 



Ho. 



f (Ha) ' 
J cont 



(2) 



where Fn a has units of erg s 1 cm 2 , Wa a has units of 
A and / c ^t is the continuum flux density at 6563 A in 
erg s _1 cm" 2 A -1 . 

In Equation [II the first term in square brackets is 
the total Ha + [N II] A6584 flux and the second term 
m square brackets is the total [N II] A6584 flux. 4 The 
measured flux densities and widths of the F658N and 
F660N filters are given by /f658Nj /f660N &nd AApgsgN, 
AAf660N, respectively. The flux densities of the contin- 
uum values, / c ^ n "' and / C ^|P, are found by interpolating 
between the F550M and i-band data points. 

Table [3] lists the measured Ha fluxes and equivalent 
widths obtained from the photometry and from direct 
measurements of the Ha line in the spectra using the 
IRAF task splot. These values agree within the er- 
rors, confirming that the appropriate imaging data can 
be used to obtain reliable Ha measurements. 

Table [3] also lists extinction corrected Ha fluxes using 
internal extinctions of A y=0.4 for Source 26 and Ay=0.5 

Galactic foregroun d extinction 
V) = 0.019, |Schlegel et 



for Source 30 (see (4.1 1. 
is also a ccounted for (E(B 



al. |1998[ ). A comparison of th e Ha measurements with" 
radio continuum observations (j Reines et al.||2008a[ ) sug- 
gest that these Ay may actually underestimate the true 
extinctions of these infant clusters since optical observa- 
tions may not probe the deepest embedded regions of the 
YMCs. Nonetheless, we use the optically derived extinc- 
tions here. 

3. PHYSICAL PROPERTIES OF THE CLUSTERS DERIVED 
FROM Ha EMISSION 

3.1. Ionizing Luminosities 

The flux and equivalent width of the Ha emission line 
can be used to infer properties of a YMC such as the 
number of ionizing photons produced by the massive 
stars, the cluster's age, and its m ass (e.g. |Alonso-Herrero 



et al. 1996 



luminosities 
minositics 



Reines et al.||2008b[ ). We calculate ionizing 
JNL yc ,_using the extinction corrected Ha lu 



from Condon 



{Table [1 
^( |1992f 



in the following equation derived 
for a 10 4 K gas: 



Lyc 



> 7.87 x 10 1 



U 



erg s~ 



(3) 



4 The [N II] line at 6548 A is a relatively weak contaminant to 
the Ho flux and the filter throughput drops at this wavelength. 
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TABLE 3 

Ha Emission and Derived Properties of the YMCs 





Observed Flux 
(erg s — 1 cm -2 ) 


Extinction Corrected Flux 
(erg s -1 cm -2 ) 


Equivalent 
Width (A) 


N Lyc 
(10 49 s" 1 ) 


Age 
(Myr) 


Mass 
(10 3 M ) 


Source 26 














Photometry 


1.1(0.1) X 10~ 12 


1.6(0.2) X 10~ 12 


1600(370) 


230(40) 


(0.1)2.9(0.8) 


65(12) 


Spectroscopy 


9.5(0.3) x 10" 13 


1.4(0.1) x 10~ 12 


1500(50) 


200(30) 


(0.1)3.0(0.1) 


62(9) 


Source 30 














Photometry 


7.4(0.9) x 10~ 14 


1.2(0.1) x 10~ 13 


3160(730) 


17(3) 


< 2.6 


4(1) 


Spectroscopy 


7.3(0.2) x 10" 14 


1.2(0.1) x 10" 13 


3260(120) 


17(3) 


< 0.4 


4(1) 



Note. — Uncertainties are shown in parentheses. 
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Fig. 5. — Starburst99 mo del e volutionary track for the equivalent 
width of Ho emission (see j ]4.1| for a description of the model). The 
measured Ha equivalent widths of Sources 26 and 30, from both 
the photometry and spectroscopy, are shown indicating their young 
ages. Uncertainties are reflected in the sizes of the shaded regions. 



Adopting a distance of 3.9 Mpc to NGC 4449, the 
production rates of ionizing photons are approximately 
2 x 10 51 and 1.7 x 10 50 s" 1 , or the equivalent of ~200 and 



17 07.5 V stars ( Vacca et al.| [l996) for Sources 26 and 
30, respectively. We note that the ionizing luminosity of 
Source 30 derived here is lower than that in lReines et al.l 
(2008al (2.7 x 10 50 s _1 ). In that paper, the ionizing lumi- 



nosity was derived from radio continuum measurements 
which probe higher extinction regions than the Ha line 
used here. 

3.2. Ages 

Ages of the YMCs are estimated from their Ha equiv- 
alent widths. The equivalent widths from both the pho- 
tometry and spectroscopy are plotted in Figure [5] with 
a Starburst99 model as a function of age. We note that 
the model equivalent widths include a nebular compo- 
nent in the continuum (as do our measurements). The 
Ha equivalent width of Source 26 suggests the cluster 
is w 3 Myr old. Source 30 has an Ha equivalent width 
consistent with the maximum value, and therefore the 
youngest ages, of the model evolutionary track. Source 



30 is certainly younger than Source 26 and is likely < 1 
Myr old. 

3.3. Masses 

Assuming mass scales with ionizing luminosity, masses 
of the clusters can be estimated by multiplying an in- 
put model mass by the ratios of the measured ionizing 
luminosities to the model ionizing luminosities at their 
res pect ive ages. Using the Starburst99 model described 

and ~ 4 x 10 3 M 
Again, the mass of 



in (4.1 



we infer masses of ~ 6 x 10 4 
for Sources 26 and 30, respectively. 
Source 30 is lower than that presented in |Reines et al 



(2008a 



(~ 6 x 10 3 M Q ) which used the radio derived ion- 
izing luminosity. Ionizing luminosities, ages, and masses 
of Sources 26 and 30, based on the Ha measurements pre- 
sented in |2.3.2[ are summarized in Table [3] It is noted 
that caution must be applied when interpreting the phys- 
ical properties of observed clusters derived from popula- 
tion synthesis models, since an underlying assumption of 
every Simple Stellar Population (SSP) model is a well- 
populated Initial Mass Function (IMF) and stochastic 
fluctuations can bias the results. 

4. MODELING THE STARBURST SPECTRAL ENERGY 
DISTRIBUTIONS 

In addition to using hydrogen recombination lines such 
as Ha, broad-band photometric observations of extra- 
galactic massive star clusters are commonly compared 
to models of SSPs in order to estimate their physical 
properties such as age, extinction and mass. Here we 
perform SED fitting to our "photom etric data using the 



Starburst99 (|Leitherer et al. 1999[) and GALEV K& 



tulla et al. 2009) evolutionary synthesis models, and in 



vestigate how well the model SEDs fit the continua in 
our observed spectra. In addition to the stellar light, 
both sets of models include nebular continuum (although 
these are computed in slightly different ways). The 
GALEV models also have the option to include a set of 
metallicity-dependen t emission lines (Anders & Fritze- 
v Alvenslebenp003 l. 

Physical properties of the clusters are estimated by 
comparing a grid of model SEDs (Starburst99 and 
GALEV separately) to measured photometric flux den- 
sities. The grid includes SEDs with ages < 50 Myr in 
steps of 0.1 Myr, and visual extinctions < 3 in incre- 
ments of 0.1 magnitudes. We apply a 30 Dor adus ex- 
tinction curve ad opted from Misselt et al. ( |1999[ Table 3) 
and Fitzpatrick ( 1985 Table 6), using the parameteriza- 
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Fig. 6. — Best-fitting Starburst99 (top) and GALEV (bottom) model SEDs for Sources 26 (left) and 30 (right) using all of the photometric 
bands (green dots) in the fitting process. Both models include stellar continuum (dashed lines) and nebular continuum (dotted lines). The 
GALEV models also include emission lines. Model parameters (masses, ages and extinctions) are shown in the upper right corners of the 
plots. The goodness-of-fit parameter, <r, is also shown for each model. A lower value of <r indicates a better fit. 

include the F170W, F336W, u, F550M, and i filters. De- 
tails about the models and the results of SED fitting are 
described below. 



(19901. Galactic fore- 



tion given by Fitzpatrick & Massa 

ground extinction tow ards NGC "4-449 is a lso accounted 
for (E(B-V) = 0.019,|Schlegel et al |1998 j) using the ex- 
tinction curve of Cardelli et al.| (fl989). Each model SED, 
of a given age and Ay, is then convolved with the total 
system throughput curves for the photometric filters to 
obtain synthetic flux densities. The best-fit model SED 
is determined by minimizing a goodness-of-fit parame- 
ter, er, equal to the standard deviation of the logarith- 
mic residuals of the observed and synthetic flux densities 
(weighted by the errors of the observed flux densities). 
A mass estimate is obtained by scaling the model mass 
by the mean logarithmic offset between the observed and 
best-fit model flux densities. 

In all, we ran our SED fitting routine four times 
for each cluster: 1) using Starburst99 and all available 
broad-band photometric data, 2) using Starburst99 and 
only photometric data with minimal line emission, 3) us- 
ing GALEV (with emission lines) and all available pho- 
tometric data, and 4) using GALEV (without emission 
lines) and only photometric data with minimal line emis- 
sion. Stellar and nebular continuum emission is included 
in the models for each case. Filters with minimal line 
e mission were identified using our Figure [2] and Table 1 



Anders & Fritze-v. Alvensleben (20031 (for the short- 



wavelength filters not covered by our spectra). These 



4.1. Starburst99: Stellar and Nebular Continuum 

We use the latest Starburst99 models (Version 5.1) 
with a me tallicity of Z = 0. 004 (as appropriate for 
NGC 4449, Lequeux et al.|1979 l, an instantaneous burst 
of 10 4 M Q with a Kroupa IMF (0.1 - 100 M ), the 
Geneva evolutionary tracks with high mass loss, and the 
Pauldrach/Hillier atmospheres. The models include stel- 
lar and nebular continuum but no emission lines. 

SED fitting is first performed using all of the avail- 
able broad-band photometric data in this study (F170W, 
F336W, u, F435W, g, F550M, r, i, F814W, and z). The 
best-fitting model SEDs for Sources 26 and 30 are shown 
at the top of Figure [6j It is clear that the model SEDs do 
a poor job fitting the spectral continua. This is because 
the spectra contain strong emission lines which are not 
accounted for in the models. The photometric bands con- 
taining significant line emission cause the model SEDs to 
be higher than the observed continua. 

Next, we perform SED fitting using only those filters 
with minimal line emission (F170W, F336W, u, F550M, 
and i). In this case, the best-fit models do a very good 
job fitting both the photometry included in the fitting 
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Fig. 7. — Same as Figure [6] except that only photometric bands with minimal line emission are included in the fitting process. Also, 
emission lines are not included in the GALEV models used here. 



process and the observed continua of the spectra 5 (see 
Figure [7| . Also, the inferred physical properties of the 
clusters (shown in Figure |7]) are consistent with those 
derived from the Ha emission (see Table p3| . 

Comparing the best-fit Starburst99 moael parameters 
from the two fitting processes (all bands versus bands 
with minimal line emission) , we find that mass estimates 
can vary significantly. For Source 30, the mass derived 
from fitting all of the photometry is a factor of 3 times too 
high due to strong lines raising the overall SED. The ex- 
tinction estimate for Source 30 is also higher by 0.7 mag 
when all of the bands are fit since the model SED needs 
to be flatter in the red to better match the F814W and 
z-band data (which contain significant line emission). 

The impact of nebular continuum on the total SEDs is 
striking in both clusters (see Figure [7]) , especially in the 
U - and /-band regions of the spectra where the Balmcr 
(3646 A) and Paschen (8207 A) discontinuities are lo- 
cated. These discontinuities, or "jumps" , in the nebular 
continuum emission are due to hydrogen and helium free- 
bound recombination processes. The nebular continuum 
actually dominates the total continuum in these wave- 
length regions for the younger YMC, Source 30. 

To further illustrate the importance of nebular contin- 
uum in young stellar populations, we plot the ratio of 

5 The flux calibration of t he blu e end of the observed spectrum 
of Source 30 is in error (see j ]2.3.1| . However, this does not affect 
our results since the photometry is being fit. 



nebular continuum to the total continuum as a function 
of wavelength for Starburst99 model SEDs (Z = 0.004) 
at various ages (Figure [8j. It is clear that, at young ages, 
nebular continuum is a significant component of the total 
continuum, even at optical wavelengths. 

4.2. GALEV: Continuum Plus Emission Lines 

GALEV models were kindly provided to us by Ralf 
Kotulla and the GALEV team. These models were run 
using the Geneva evolutionary tracks with a minimum 
age and time resolution of 0.1 Myr, unlike the models 
available on the GALEV website which use the Padova 
isochrones and have a minimum age and time step of 4 
Myr. The models also use a "fraction of visible mass" 
equal to 1 rather than 0.5, which is used for the stan- 
dard models available on the web (Ralf Kotulla, private 
communication) . The models were run with a metallicity 
of Z = 0.004 and a Kroupa IMF (0.1 - 100 M ). 

These models are very similar to the Starburst99 mod- 
els described above, with the notable exception that the 
GALEV models include emission lines from nebular gas 



(Anders & Fritze-v. Alvensleben 2003 Kotulla et al 



2009p . The flux of the hydrogen lines are computed using 



atomic physics and the production rate of ionizing pho- 
tons, whereas non-hydrogen line strengths are computed 
using metallicity-dependent line ratios relative to H/3. In 
the case of low-metallicity gas (including Z = 0.004), 
observed line ratios are taken from the large sample of 
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Fig. 8. — Ratio of the nebular continuum to the t otal continuum 
as a function of wavelength for a Starburst99 SSP (j ]4,l| at various 
ages. At young ages, nebular continuum is a significant component 
of the total continuum at optical wavelengths, especially shortward 
of the Balmer (3646 A) and Paschen jumps (8207 A) located in the 
U- and /-band wavelength regions. The impact of emission lines is 
not included here. 



giant H II region s and Blue Compact Dwar f (BC D) galax- 
ies presented in Izotov et al. ( 1994 1997[ ) and Izotov & 
Thuan| ( [l998l ). 

We perform SED fitting with the GALEV models in 
the same way as the Starburst99 models. First we fit 
all of the photometric filters and include emission lines 
in the models (the option to include lines is not avail- 
able for Starburst99). The best- fitting model SEDs for 
Sources 26 and 30 are shown in Figure [6] It is clear that 
the GALEV models including emission lines do a much 
better job matching the continua in the observed spectra 
than the Starburst99 models do when we fit all of the 
photometry. However, the GALEV model SEDs shown 
in Figure [6] are still a bit too high in the red to near- 
infrared and too low in the ultraviolet. This is likely 
due to mismatched line intensities between the mod- 
els and the observed spectra. We note that currently 
the GALEV models are missing the higher level Paschen 
lines, although they do include Paa through Pae. Over- 
all, however, the GALEV models including emission lines 
perform much better than the Starburst99 models when 
we fit all of the photometry including bands with signif- 
icant line contamination. 

We also only fit the photometry with minimal line 
emission with GALEV SEDs as we did with the Star- 
burst99 models. Here we exclude emission lines in the 
models so they only contain components from the stellar 
and nebular continua. The best-fitting SEDs are shown 
in Figure [7J These models perform equally well as the 
Starburst99 models fit only to the photometry with min- 
imal line emission (also shown in Figure [7]). However, 
it is curious that the best-fit model parameters from the 
two families of models are as different as they are. The 
differences may be due to the ways in which the nebu- 
lar continuum is calculated. The best-fit GALEV SEDs 
shown in Figure [7] also do a better job fitting the con- 
tinua of the observed spectra than the GALEV models 
(including lines) fit to all of the photometry (Figure |6| . 
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Fig. 9.— Spectrum of Source 26 multiplied by the HST/ACS 
WFC F814W total system throughput curve. The best-fitting Star- 
burst99 model SED from FigurelTlis shown as a solid blue line. The 
stellar and nebular continua are shown as dashed and dotted lines, 
respectively. The relative contributions of the stellar continuum, 
nebular continuum and emission lines to the total flux of this filter 
are approximately 40%, 40% and 20%. 

5. THE IMPORTANCE OF NEBULAR EMISSION IN 
BROAD-BAND PHOTOMETRY 

A number of authors have discussed the impact of neb- 
ular emission on broad-band photometry of you ng stellar 
populations based on synthesis modeling (e.g. |Krueger 
Leitherer fc Hcckman |1995| fZackriss on et 



et al 



al 



al 



199o| , 

2 001I |Anders fc Fritze-v. Alvensleben||2003[ |Molla et 



2009) 



Our spectrophotomctic observations of infant 
massive star clusters clearly support the conclusion that 
nebular emission, both in lines and continuum, signifi- 
cantly affects broad-band photometry of these objects. 

5.1. Contributions to Broad-band Fluxes 

In an effort to provide a benchmark for estimating 
the impact of nebular emission on photometric observa- 
tions of young massive stellar populations, we compute 
the relative contributions of the stellar continuum, neb- 
ular continuum, and emission lines to the total flux of 
Source 26 through various HST filter/instrument com- 
binations (including WFC3) and the SDSS griz filters. 6 
The contribution of nebular emission in the HST/ACS 
WFC F814W filter is of particular interest sinc e, in their 
stu dy of ra dio-detected clusters in NGC 4449, Reines et 
al. ( 2008a ) found an excess in this band relative to model 
SEDs for a pure stellar continuum. 

The spectrum of Source 26 multiplied by the HST/ACS 
WFC F814W total system throughput curve is shown in 
Figure 9 The best-fitting Starburst99 model SED (from 
Figure 7 1 is also shown. In addition to stellar and nebu- 
lar continuum, the spectrum contains several prominent 
emission lines in this band. The strongest are [S III] 
A9532 and [S III] A9069. A number of other lines, in- 
cluding the Paschen series, are detected and labeled in 
Figure [9} 

We determine the relative contributions of stellar con- 
tinuum, nebular continuum and emission lines to the to- 

6 We do not use the spectrum of Source 30 since it has a low 
signa l-to- noisc ratio and the flux calibration is uncertain in the blue 
(see ^TJand p.3.1\ . 
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TABLE 4 

Relative Contributions to Broad-band Fluxes of Source 26 



Filter Instrument Stellar Nebular Emission 

Continuum Continuum Lines 



HST Medium and Wide Filters in the /-band Region 



F845M 


WFC3/UVIS 


0.49 


0.41 


0.09 


F814W 


WFC3/UVIS 


0.39 


0.38 


0.23 


F814W 


ACS/WFC 


0.41 


0.39 


0.20 


F814W 


ACS/HRC 


0.38 


0.36 


0.25 


F814W 


WFPC2 


0.41 


0.40 


0.20 


F791W 


WFPC2 


0.46 


0.45 


0.10 


F775W 


WFC3/UVIS 


0.46 


0.45 


0.09 


F775W 


ACS/WFC 


0.46 


0.46 


0.08 


F775W 


ACS/HRC 


0.46 


0.46 


0.08 


F763M 


WFC3/UVIS 


0.47 


0.50 


0.03 


HST Medium and Wide Filters in the V-band Region 


F555W 


WFC3/UVIS 


0.31 


0.12 


0.56 


F555W 


ACS/WFC 


0.27 


0.11 


0.63 


F555W 


ACS/HRC 


0.26 


0.11 


0.63 


F555W 


WFPC2 


0.32 


0.13 


0.54 


F550M 


ACS/WFC 


0.70 


0.31 


0.00 


F550M 


ACS/HRC 


0.70 


0.31 


0.00 


F547M 


WFC3/UVIS 


0.71 


0.30 


0.00 


F547M 


WFPC2 


0.68 


0.29 


0.03 



SDSS Filters 



SDSS 0.32 0.10 0.58 

SDSS 0.27 0.16 0.57 

SDSS 0.46 0.46 0.07 

SDSS 0.34 0.26 0.40 



tal observed flux in the F814W filter as follows. The 
contributions from the stellar and nebular continua are 
estimated by multiplying the models for these compo- 
nents with the filter throughput curve and integrating 
over the bandpass. The total flux from emission lines 
in this filter is determined by multiplying the observed 
spectrum with the F814W filter throughput curve and 
measuring the flux of each line in the resulting spectrum. 
The relative contributions of each component arc found 
by dividing by the total observed flux in the F814W fil- 
ter (obtained by convolving the observed spectrum with 
the throughput curve). We find that the relative con- 
tributions of the stellar continuum, nebular continuum 
and emission lines to the total observed flux of Source 26 
through the F814W filter are approximately 40%, 40%, 
and 20%, respectively. Clearly, nebular emission con- 
tributes a significant fraction (~ 60%) of the total flux 
in the F814 W filter and should not be neglected. The 
combination of nebular continuum and emission lines can 
acc ount for the F814W /-band excess found by |Reines et| 



al. (2008a). This result is consistent with their predic- 
tion matTtne origin of the excess should not affect clusters 
older than ~ 5 Myr. 

Relative contributions to the total flux of Source 26 
through other HST filter/instrument combinations are 
listed in Table [3] Throughput curves were obtained us- 
ing SYNPHOT, the synthetic photometry package dis- 
tributed by the Space Telescope Science Data Analysis 
System (STSDAS). The total flux and the flux of the 
stellar and nebular continua are found in the same way 
as described above. The line contributions, however, are 



simply estimated by subtracting the continuum contribu- 
tions from the total observed flux. Plots of the observed 
spectrum of Source 26 multiplied by the HST WFC3 fil- 
ters listed in Table [4] are shown in Figure [lO] We find 
that the relative contribution of the nebular gas (contin- 
uum plus lines) to the total flux of Source 26 in any of the 
HST /-band filters listed in Table g] is at least 50%. In 
the U-band, roughly 70% of the F555W flux comes from 
ionized gas. Younger clusters will have an even larger 
contribution from ionized gas. Relative contributions to 
the SDSS fluxes of Source 26 are also listed in Table |4] 
using the throughput curves available online. 7 

We emphasize that the quantitative results listed in 
Tableware specific to Source 26. Relative contributions 
will vary from source to source, depending on (at a mini- 
mum) age and metallicity (e.g. our Figure [8] and Anders 
& Fritze-v. Alvensleben 2003). However, it is evident 



that the impact of nebular emission on broad-band fluxes 
can be quite large and should not be readily dismissed. 

5.2. The Impact on Magnitudes, Colors and Derived 
Properties 

The ionized gas associated with a young stellar popu- 
lation can provide a significant cont ribution to the total 
observed flux in a given filter ([5.1|. The nebular emis- 
sion, therefore, can have large impact on the physical 
properties derived from magnitudes and broad-band col- 



Using the GALEV models described in Section |4.2[ we 
have computed the magnitude difference as a function 

7 http:/ /www. sdss.org/dr6/instruments/imager/ 
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Fig. 10. — Spectrum of Source 26 multiplied by the HST/WFC3 filters in the V- and /-band wavelength regions listed in Table [4] The 
best-fitting Starburst99 model SED from Figure u\ is shown as a solid line. Individual contributions from the stellar and nebular continua 
are shown as dashed and dotted lines, respectively. The impact of emission lines on the flux through various filters is given in Table [4] The 
F555W and F814W fluxes are most affected by line emission. 



of age between a pure stellar population and one of the 
same mass but including ionized gas emission (contin- 
uum plus lines) for the HST/WFC3 F814W (broad I), 
F555W (wide V), and F547M (narrow V) filters. The 
results are shown in Figure [IT] For a 3 Myr old clus- 
ter, including the gaseous emission results in magnitudes 
that are brighter by 1.0 (F814W), 1.5 (F555W), and 0.4 
(F547M) mag. Neglecting nebular emission and assum- 
ing all of the flux comes from the stars would therefore 
lead to overestimates of the stellar mass by factors of 2.5 
(F814W), 4.0 (F555W), and 1.4 (F547M). From Figure 
we can see that significant contributions from both 
nebular continuum and lines are present in the F814W 



filter, whereas the F555W flux is most affected by strong 
lines (e.g. H/3 and [O III]). The F547M flux is affected by 
the nebular continuum but does not contain strong line 
emission. 

We obtain similar results from modeling the SED of 
Source 26 and determining the relative contributions of 
the stellar continuum, nebular continuum and emission 



lines to the total flux through various filters ({ 5.1 1. Mass 
(over)estimates of this YMC (age » 3 Myr, 0.004), 
assuming all of the flux in a given filter is due only to 
stars, can be obtained by taking the inverse of the relative 
contribution of the stellar continuum to the total flux 
(the third column in Table [4]). For example, neglecting 
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Fig. 11. — Left: Magnitude difference as a function of age between a pure stellar populatio n an d one of the same mass but including 
ionized gas emission (continuum and lines). The GALEV model SEDs described in Section |4.2| were convolved with the HST/MVFC3 
F814W (broad /), F555W (wide V), and F547M ( narrow V) niters to p roduce this plot. Rigfti:Color-color diagram using the HST/ ACS 
WFC F435W, F550M and F814W filters and the data of |Reines et al | |2008a) . Sources 26 and 30 are shown as large purple and orange 
dots, respectively. Model evolutionary tracks are shown for a Starburstyy model with stellar continuum only (yellow), a Starburstyy model 
with stellar and nebular continuum (blue), and a GALEV model with stellar continuum, nebular continuum, and emission lines (red). The 
tracks extend to an age of 50 Myr. An age of 1 Myr is indicated by a star and the dots along the tracks indicate ages of 2, 3, 4, 5, 7, 10, 
25 and 40 Myr. Observed and model colors are for the HST filters with VEGAMAG zeropoints. 



the effects of nebular emission in the ACS/ WFC F555W 
filter would result in a mass 3.7 x larger than the mass 
obtained by properly accounting for the gas contribution 
in this filter. 

Since the magnitudes of young stellar populations can 
be considerably brightened by the presence of ionized 
gas, colors can als o be significantly affe cted. For ex- 



ample, the work of Reines et al. (2008a I found an ap- 
parently anomalous color-color diagram in which the 
radio-detected YMCs in NGC 4449 with the bluest 
HST/ ACS WFC [F435W-F550M] colors had the red- 
dest [F550M-F814W] colors, contradictory to the Star- 
burst99 model evolutionary track for a pure stellar con- 
tinuum. We have reproduced this color-color diagram 
but have now included a Starburst99 model evolution- 
ary track including nebular continuum, and a GALEV 
track i nclu ding nebular continuum plus emission lines 
(Figure 11). Only the GALEV model is able to repro- 
duce the trend seen in the data, signifying that emission 
lines have a large effect on these colors at young ages (< 5 
Myr). Although the F550M filter is essentially line-free, 
the F814W and F435W filters contain a number of emis- 
sion lines. The forbidden [S III] lines at 9532 and 9069 A 
are the most prominent in the F814W filter, but many 
others are also present (see Figure 9l). The F435W filter 
contains line emission from H/y, IM, He+[N III] A3970, 
HC+[He I] A3 889, [N III] A3869, and [O II] A3727 (se e our 
Figure [2] and |Anders fc Fritze-v. Alvensleben |2003 1 . 

It is worth noting that the [F435W-F550MJ versus 
[F550M— F814W] color-color diagram may prove to be 
a useful diagnostic for identifying YMCs with ages < 5 
Myr. It is also interesting that the extinction vector in 
Figure [TT] is perpendicular to the GALEV model track 
for ages % 5 Myr. In other words, there is no degeneracy 
in this color space between age and extinction at young 
ages. Therefore, rough estimates of both age and extinc- 
tion could be obtained from this color-color plot. How- 



ever, we emphasize that the line strengths of any given 
source may vary from the ones input into the GALEV 
models which could potentially lead to large uncertain- 
ties in the derived properties. 

As a check, we look at Sources 26 and 30 (indicated in 
Figure [TT] by purple and orange dots respectively) . The 
ages inferred from their colors are consistent with those 
obtained from their Ha equivalent widths and from SED 
fitting; approximately 3 Myr for Source 26 and 1-2 Myr 
for Source 30. The visual extinction of Source 26 also 
appears to be consistent with our previous estimates of 
^0.4 magni tud es. The extinction of Source 30 inferred 
from Figure 11 (« 1 mag), on the other hand, is half a 
magnitude larger than that obtained from SED fitting to 
the photometry with minimal line emission (Figure M). 8 
This difference is likely the result of different line inten- 
sities between the models and the data. The GALEV 
models, however, at least qualitatively explain these col- 
ors of very young star clusters, which neither of the Star- 
burst99 models can do. 

6. CONCLUSIONS 

This paper examines the importance of nebular con- 
tinuum and line emission in observations of young mas- 
sive stellar populations. We have obtained spectroscopy 
(A ~ 3800 - 9800 A) of two young (< 3 Myr) star clus- 
ters in the nearby starburst galaxy NGC 4449 and sup- 
plement these data with archival HST and SDSS pho- 
tometry of the clusters. We compare our data to the 
Starburst99 and GALEV evolutionary synthesis models 
to estimate the physical properties of the clusters and 
determine the impact of ionized gas emission (line and 
continuum) on broad-band photometry. Our main re- 
sults are summarized below: 

8 It is not clear that optically derived Av 's are accurate measures 
of the extinction at such young ages, since a significant fraction of 
the cluster may be heavily embedded and visibly obscured. 
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The combination of nebular continuum and line 
emission can account for the F814W /-band ex- 



cess found by Reines et al. (2008a I in their study 
of radio-detected YMCs in NGC 4449. This result 
is consistent with their prediction that the origin 
of the excess should not affect clusters older than 
~ 5 Myr. 

2. At very young ages (< 3 Myr), nebular contin- 
uum emission from ionized gas can rival the stellar 
luminosity of YMCs at optical wavelengths. The 
relative contribution from the nebular continuum 
emission is largest in the U- and /-bands, where 
the Balmer (3646 A) and Paschen jumps (8207 A) 
are located. 

3. Nebular line emission is significant in many com- 
monly used broad-band HST filters including the 
F814W /-band, the F555W U-band and the 
F435W Z?-band. Emission lines (in addition to neb- 
ular continuum) are required to explain the dis- 
tribution of radio-detected YMCs in NGC 4449 
in the [F435W-F550M] versus [F55 0M-F814W] 
color-co lor diagram first presented in |Reines et al | 
( |2008aD . 

4. SED fitting to broad-band photometry of YMCs 
is most successful when the photometric bands do 
not contain significant line emission. In this case, 
evolutionary synthesis models only containing con- 
tinuum emission from stars and ionized gas are ade- 
quate. If the photometric bands include significant 



line emission, however, models including emission 
lines perform much better than those that do not. 

5. The impact of ionized gas emission on broad-band 
fluxes, magnitudes and colors can be large and 
should not be readily dismissed. We urge caution 
when comparing observations of young (< 5 Myr) 
clusters to synthesis models since nebular contin- 
uum and line emission can significantly affect in- 
ferred properties such as ages, masses and extinc- 
tions. 
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